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THEORY  AND  APPLICATIONS  OF  ELECTRON  TRANSFERS  AT  ELECTRODES  AND  IN  SOLU- 
TION 

R.  A.  MARCUS 

Department  of  Chemistry,  University  of  Illinois,  Urbana,  Illinois  61801 


ABSTRACT 


The  theory  of  simple  electron  transfers  at  electrodes  and  in  solution 
is  reviewed,  and  various  thermal  fluctuations  of  coordinates  leading  to 
electron  transfer  are  described.  A simplified  derivation  is  given  of  the 
free  energy  of  such  fluctuations.  Implications  of  the  theoretical  equa- 
tions for  experiment  are  describe^,'  some  of  the  relevant  results  having 
been  summarized  recently  in  Dahlem  Konferenz  Phys.  Chem.  Sc i . Res.  Rept., 

1(  1975)  477.  They  inciude  relations  between  rate  constants  of  cross- 
reactions and  se I f-exchange  reactions,  between  rates  of  reactions  at  elec- 
trodes and  those  in  solution,  nonspecific  solvent  effects,  chemilumines- 
cence, and  other  properties.  Approximate  equations  of  the  BEBO  type  are 
also  given  for  Tafel  slopes,  Bronsted  slopes  and  rates  of  cross-reactions, 
for  systems  involving  rupture  and  formation  of  bonds. 
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INTRODUCTION 


The  theory  of  reactions  at  electrodes  has  much  in  common  with  that 
of  reactions  in  solution.  The  electrode  behaves  as  one  large  reactant, 
one  with  special  properties:  it  has  numerous  electronic  energy  levels  and 

the  energy  of  those  levels  is  controllable  by  the  electrode  potential. 

Just  as  one  reactant  in  solution  may  bind  the  other,  the  electrode  may 
adsorb  the  other  reactant. 

In  the  present  paper  we  shall  be  concerned  principally  with  simple 
electron  transfer  reactions  at  electrodes  and  in  solution,  reactions 
which  involve  no  rupture  or  formation  of  chemical  bonds.  Elsewhere  [I] 
we  have  discussed  reactions  involving  rupture  of  a chemical  bond  at  an 
electrode  M,  e.g.,  Eqs.  [I] -[2],  and  for  brevity  omit  a discussion  of  this 
topic  here. 
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A simple  electron  transfer  reaction  at  an  electrode  can  be  written 
as 

A,(ox)  + M(ne)  - A ,(  red)  + M , [3] 

while  a simple  electron  transfer  between  species  A,  and  A2  in  solution  is 
represented  as 

A i(  Ox)  + Aj(  red)  -*  A |(  red)  + A2(  ox)  , [4] 

where  ( ox)  and  (red)  denote  the  oxidized  and  reduced  forms  of  the  two 
chemical  species  At  and  A2.  The  ionic  charges  on  these  reactants  are 
usually  about  one  to  three,  though  sometimes  as  high  as  four  and  some- 
times as  low  as  zero.  The  solvation  energies  are  therefore  very  large 
and  their  fate  during  the  course  of  the  reaction  must  be  analyzed  care- 
fully. 

POTENTIAL  ENERGY  CURVES,  SURFACES,  FLUCTUATIONS,  AND  RATES 

It  is  useful  to  examine  first  the  vibrational  motion  within  the  re- 
actant in  Eq.  [3],  considering  initially  the  case  of  one  vibrational  co- 
ordinate q.  (The  same  plot  suffices  for  Eq.  [4]  also.)  A plot  of  the 
potential  energy  U versus  q for  the  system  on  the  left  hand  side  of  Eq. 

[3]  is  labelled  R in  Fig.  I,  and  a plot  of  U versus  q for  the  system  on 
the  right  hand  side  of  Eq . [3]  is  labelled  P [2,3]. 

There  are  many  electronic  states  of  the  metal--a  continuum  of  them. 

Thus,  Fig.  1 for  the  reaction  in  Eq.  [3]  should  consist  of  many  parallel 
curves,  vertically  displaced  from  each  other,  one  for  each  electronic 
state.  However,  one  can  show  [3,4]  that  the  electrons  donated  from  the 
electrode  M to  A( ox)  in  Eq.  [3]  come  from  levels  near  the  Fermi  level, 
and  most  of  them  donated  from  A( red)  to  M in  the  reverse  reaction  go  into 
levels  near  the  Fermi  level.  Thus,  it  suffices  to  confine  our  attention 
for  the  present  purposes  to  the  two  curves  in  Fig.  I and  take  the  energy 
level  of  the  electron  in  the  metal  as  the  Fermi  level.  The  following 
analysis  then  applies  equally  well  to  reactions  [3]  and  [4], 

Several  facts  are  noted: 

(1)  The  minima  of  the  two  curves  occur  at  different  values  of  q,  re- 
flecting the  fact  that  the  equilibrium  bond  length  in  A is 

ox 

• different  (usually  shorter  in  the  case  of  transition  metal  ions 

and  a metal-ligand  bond)  from  that  in  A 

(2)  The  relative  height  of  the  two  minima  LUCI  depends  on  the  elec- 
trostatic potential,  the  P curve  being  lowerf;d  vertically 

_J 


POTENTIAL  ENERGY 


D 
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Fig.  1.  A plot  of  the  potential  energy  U of  the  system  consisting  of 
reactants  plus  solvent  (R),  along  some  coordinate  q,  and  of  the  system 
consisting  of  products  plus  solvent  ( P) , holding  all  other  coordinates 
fixed,  for  reaction  [3]  or  [4]. 

relative  to  the  R curve  by  making  the  electrode  more  negative, 
i.e.,  by  decreasing  e(cp^^p^)  where  the  cp's  represent  the  po- 
tentials of  the  met3l  and  solution. 

(3)  When  the  difference  of  equilibrium  bond  lengths  £q°  is  less,  fo 
a given  AU°,  the  barrier  height  AU*  is  less. 

(4)  When  the  P curve  is  lowered  relative  to  the  R one,  the  barrier 
height  All'”",  corresponding  to  the  intersection  of  the  R and  P 
curves,  is  reduced. 

The  two  curves  have  further  properties:  When  the  reactant  is  far 

from  the  electrode  the  appropriate  curves  in  Fin.  I merely  cross  in  the 
intersection  region  (dotted  lines  there).  When  the  reactant  is  close  to 
the  electrode,  the  electronic  interaction  of  the  orbitals  of  the  reactant 
and  the  electrode  perturbs  the  dotted  line  curves,  particularly  near  the 
point  of  intersection.  At  that  point  the  unperturbed  R and  P electronic 
quantum  states  are  degenerate,  and  the  degeneracy  is  broken  by  the  inter- 
action. The  new  curves  are  the  solid  curves.  The  energy  of  the  maximum 
of  the  lower  solid  curve  is  less  than  that  at  the  intersection  by  an 
electronic  interaction  energy  S|2, 


the  "resonance  energy." 
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When  G i 2 is  zero  (e.g.,  large  separation  distance  between  electrode 
and  reactant)  a system  having  a potential  energy  curve  R will  undergo 
thermal  fluctuations  of  the  vibrational  coordinate  q,  but  not  transfer 
to  the  P curve  even  on  passing  through  the  intersection  region.  When  e , 2 
is  nonzero,  there  is  a coupling  of  the  electronic  orbitals  of  the  react- 
ant and  the  electrode,  and  there  is  a nonzero  proDability  x of  going  from 
the  R curve  to  the  P curve  when  the  system  passes  through  the  intersec- 
tion region.  When  e ) 2 is  small  the  reaction  has  a small  x and  the  reac- 
tion is  termed  nonad i abat ic.  When  x is  of  the  order  of  unity  the  reaction 
is  adiabatic.  The  electronic  factor  x can  be  calculated  by  the  Landau- 
Zener-Stueckel berg  formula  [5],  or  by  a more  numerical  quantum  mechanical 
computation  [6],  if  e12  and  the  properties  of  the  R and  P curves  are  known. 

One  notes  from  Fig.  i that  a fluctuation  of  coordinate  q is  needed 
for  the  electron  transfer  to  occur:  Not  only  must  the  reactant -e I ect rode 

distance,  or  in  Eq.  [4]  the  reactant-reactant  distance,  become  sufficiently 
small  for  electronic  orbital  overlap  to  occur  but  also  the  transfer  does 
not  occur  at  the  equilibrium  value  of  q for  the  R-curve.  Only  by  absorp- 
tion of  light  could  the  system  move  from  the  R curve  to  the  P curve  at 
this  q.  Otherwise  the  vibrationa I momentum  p,  given  by  „'2m(  E-Uj,  would 
not  be  conserved  during  the  electronic  transition  and  thus  the  Franck- 
Condon  principle  would  be  violated.  However,  at  the  crossing  point  this 
momentum  p is  substantially  conserved, so  the  transfer  from  the  R curve  to 
the  P curve  can  occur. 

A simple  calculation  provides  an  estimate  of  the  barrier  height  AU'"' 
in  Fig.  I:  It  will  be  assumed  that  the  resonance  energy  Gi?  is  small 

enough  that  the  height  of  the  crossing-point  is  approximately  that  of  the 
maximum  of  the  lowest  solid  curve  in  Fig.  I.  We  let  the  potential  energy 
of  the  R and  P curves  be  approximated  by 

Ur(q)  = ^ k(  q -q  0 r ) 2 [5] 

UP(q)  ",-k(q-q0V  + AU°  , l6i 

where  qa'  end  q0P  are  the  equilibrium  oorid  lengths  for  the  A(ox)  and 
A( red) , respectively.  AU°  is  negative  in  Fig.  I.  and  is  a linear  func- 
tion of  ne(cp  -<D  ) in  the  case  of  reaction  f 3 ; At  the  intersection  of 

m s 

the  dotted  curves  we  have 

ur  = Up  , [I] 

an  equation  which  can  be  solved  for  the  value  of  q,  q*  at  the  intersec- 


tion.  Introducing  this  q^  into  Eq.  [bj  one  obtains  ultimately 

AU*  = ( X | / 4 ) [ I + (AU°/\.)]  2 , [8] 

where 

X.  = j k(  Aq  °) 2 , Aq 0 = qoP  - q°r  . [9] 

\./4  is  actually  the  barrier  height  when  the  AU°  in  Fig.  I is  zero. 

The  rate  constant  k^  would  be  given  by  the  collision  frequency  Z per 
unit  area  of  electrode  per  unit  time  in  the  electrode  case  and  per  unit 
concentration  per  unit  time  in  the  homogeneous  case,  multiplied  by  the 
Boltzmann  factor  exp( -AU*/RT]  , by  a Boltzmann  factor  exp(  -wr/RT)  for 
bringing  the  reactant  to  the  electrode,  and  by  h: 

k(_  = k Z exp[-(  wr+AU*)/RT]  . [10] 

Here,  wr  is  the  work  term,  if  any,  required  to  bring  the  reactants  to- 
gether to  some  separation  distance  R.  When  such  work  terms  occur,  the 

AU°  in  Eq.  [8]  is  replaced  by  All0,  the  AU  at  that  separation  distance  R. 

K 


Thus,  Eq.  [ 8]  is  replaced  by 

AU*  = (X,/4)[l  + (AUr°/X;)]  2 

[M] 

AUp°  = AU°  + wp  - wr  , 

[12] 

where  AU°  is  the  AU  at  infinite  separation  and  wP  is  the  work  to  bring 
the  products  together  to  the  separation  distance  R. 

The  above  calculation  of  AU'*  is  classical,  and  indeed  classical 
mechanics  is  commonly  used  nowadays  to  treat  reactive  collisions,  by  cal- 
culating the  trajectories  of  all  the  atoms  of  a reacting  molecule  or  pairs 
of  molecules  during  a collision.  However,  for  some  problems  a quantum 
mechanical  treatment  is  needed,  for  example  for  treating  a protonic  vibra- 
tion in  Eq.  [I],  A quantum  treatment  is  given  in  Refs.  [7,8],  other  ap- 
proximations being  introduced  as  well. 

The  results  in  Eqs.  [10] -[12]  were  derived  for  the  case  of  one  vi- 
bration. The  derivation  can  be  extended  to  all  the  vibrations  of  the  re- 

actant] s ) in  Eq.  [3]  or  [4],  The  potential  energy  curves  of  Fig.  I are 
replaced  by  potential  energy  surfaces,  plotted  as  a function  of  all  the 
q's  in  the  system  rather  than  just  one.  It  there  are  N q's.  the  inter- 
section which  is  a point  in  Fig.  I becomes  a h,  per  surface  of  ‘J-l  dimen- 
sions i n the  fj-d  i mens  i ona  I case.  The  coordinate  q in  Fig.  I then  repre- 
sents some  path  in  the  It-d  imens  > ona  I q-space  and  tin  R-  and  F-curves  arc- 
profiles  o.  the  ictual  potential  energy  suriaces  plotted  along  that  pot 


0 


1 

In  effect,  becomes  a sum  of  terms  of  the  type  in  Eq.  [9],  summed  Over 
al 1 vibrations,  i . e. , 

Xj  = j Z.k.(Aq.0)2  , [13] 

and  k.  is  related  to  the  force  constants  of  a bond  k.r  and  k.*3  in  the 
i i i 

Oxidized  and  reduced  forms  [3],  A somewhat  more  sophisticated  treatment 
of  \.  is  given  in  Ref.  [3],  leading  to  k.  being  a certain  average  of  k.f 

and  kj^. 

There  remain  the  solvent  fluctuations  outside  of  the  inner  coordina- 
tion shell  of  the  reactant  in  Eq.  [3]  or  reactants  in  Eq.  [4],  Here,  the 
potential  energy  functions  do  not  depend  on  the  solvent  coordinates  (ori- 
entations, translations)  in  the  simple  quadratic  fashion  in  Eqs.  [5]  and 
[6]  and  of  Fig.  I.  The  treatment  of  the  solvent  coordinates  is  corres- 
pondingly more  complicated.  However,  one  feature  is  immediately  clear: 

Just  as  a thermal  fluctuation  of  vibrational  coordinates  was  needed  to 
reach  the  intersection  region  in  Fig.  I,  a suitable  thermal  fluctuation 
of  solvent  orientation  coordinates  or  reactant's  vibrations  also  permits 
the  system  to  reach  the  N-l  dimensional  hypersurface  (the  intersection 
region).  A statistical  mechanical  treatment  of  the  free  energy  associ- 
ated with  these  fluctuations  is  given  in  Ref.  [3],  Dielectric  continuum 
theory  also  permits  an  estimate  of  the  latter  to  be  made,  and  was  first 
given  for  the  homogeneous  reaction  case  in  1956  [9],  and  in  ONR  Technical 
Report  Mo.  12  (presented  e'sewhere  in  this  volume)  in  1957  for  the  elec- 
trode reaction,  Eq.  [3],  i 

The  intersection  region  in  Fig.  1,  now  for  a generalized  coordinate 
q and  a many-dimensional  configuration  space,  is  reached  by  a thermal 
fluctuation  of  coordinates  q from  their  most  probable  values.  Such  a fluc- 
tuation produces  a corresponding  thermal  fluctuation  in  the  solvent  di- 
electric polarization  function.  In  the  intersection  region  this  polari- 
zation is  appropriate  neither  to  the  charges  of  the  reactants  nor  to  those 
of  the  products  but  rather  to  some  compromi  se.  wh  i ch  depends  on  All0  or  more 
generally  on  a free  energy  of  reaction  AG0'  or  in  the  electrode  case  on  a 
half-cell  potential  E minus  the  standard  potential  E0',  E — E 0 1 . This  di- 
electric polarization  was  termed  nonequilibrium  polarization  [9]  and  in- 
deed the  original  dielectric  continuum  ( nonstat i st ica I mechanical)  deri- 
vation  [ 9 , 10]  was  concerned  with  fluctuations  in  solvent  dielectric  pol- 
a r i z a t ion. 

I 


It  is  instructive  to  give  a somewhat  simpler  derivation  of  the  free 


energy  change  needed  to  reach  the  intersection  region  by  these  fluctua- 
tions of  solvent  dielectric  polarization.  This  derivation  is  given  in 
the  next  section. 
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SIMPLIFIED  DERIVATION  OF  NONEQUILIBRIUM  POLARIZATION  EXPRESSION 

We  consider  the  homogeneous  reaction  system  Eq , [4]  first  and  then 
indicate  the  modification  for  the  electrode  case  Eq.  [3].  We  denote  the 
charges  by  e.  and  radii  by  a.  for  reactants  I and  2 (i  = '.2),  and  add  a 
superscript  p to  the  e.  to  denote  charges  of  the  products.  The  static 
dielectric  constant  of  the  solvent  medium  is  denoted  by  and  the  square 
of  the  refractive  index  (the  "optical  dielectric  constant")  by  D . The 
separation  distance  is  denoted  by  R. 

A nonequilibrium  dielectric  polarization  of  the  medium  can  be  pro- 
duced in  a reversible  manner  by  a two-step  charging  process.  Since  each 
step  is  reversible,  the  free  energy  of  formation  of  this  nonequilibrium 
system,  i . e. , the  free  energy  of  this  polarization  fluctuation,  can  be 
calculated  in  a relatively  straightforward  manner.  The  two-step  charging 
process  is  the  following,  at  a given  separation  distance  R. 

(1)  Change  the  charge  of  each  reactant  i from  e.  to  e . 1 , e.‘  being 
so  chosen  to  produce  the  desired  or  ientai.  iona  1 -vi  brat  iona  I di- 
electric polarization. 

(2)  Change  the  charge  of  each  particle  i back  from  e. ' to  e.,  hold- 
ing the  above  or ientat iona I -vibrat iona I dielectric  polarization 
fixed. 

The  details  of  the  calculation  are  as  follows,  where  the  electrostatic 
potential  in  the  solvent  medium  at  any  point  r is  denoted  by  i(r). 


Step  I 

The  val„e  of  e.  and  y(  r)  at  any  stage  v of  the  charging  process  are 
denoted  by  e.v  and  ')/v(r),  respectively.  They  are  given  by 


D r 
s 


i 


[14] 


e.v  = e.  + v( e . 1 -e.)  [lb] 

i i vv  i i ' ’ 

where  r.  is  the  distance  from  the  field  point  r to  the  center  of  ion  i. 

V varies  from  0 at  the  beginning  of  the  charging  process  to  I at  the  end. 
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and  e.v  can  be  written,  thereby,  as  in  Eq.  [15], 

The  potential  at  the  surface  of  ion  I due  to  the  medium  and  to  ion  2 
is  obtained  by  replacing  r,  by  a,  inEq.  [14],  § ,v,  the  potential  there 
minus  the  seif-potential,  is  obtained  by  subtracting  e,v/a.  from  [14]: 


$ 


e„v  e v 1 

§^7  + rHf  - '> 

s ^ s 


[16] 


The  average  of  J ,v  over  the  surface  of  ion  I is  denoted  by  $ ,v  and  is 
found  to  be 


$ ■ ' 


e v _ v i 

“2 — 4 t-L-1  _L  _ 1 ) 
D R aT'o  ' 
s s 


[17] 


The  average  leading  from  Eq.  [16]  to  [17]  will  be  recognized  as  the  well- 
known  electrostatic  result  [II]  that  the  average  value  of  a 1 / r 2 poten- 
tial from  a uniform  distribution  over  a sphere  is  l/R. 

When  |,v  is  multiplied  by  an  increment  of  charge  d(e(vj,  i.e.,  (e,1- 
e,)dv,  and  integrated  over  v from  0 to  1 and  when  the  same  integration  is 
performed  for  ion  2,  and  both  terms  summed  we  obtain  the  work  term  Wj 
required  in  charging  step  I; 


Wj  = j 5 i “(  e | 1 -e  ,)  dv; 


i2v(  e2 1 -e2)dv  . 


[18] 


o o 

Eqs.  £17] — [ 18]  yield  immediately 


= + [e  ,6e  , + ^ Ae  ,) 2 j ( rj- - l)j-  + [e2Ae2  + ~{Lez)  2](:f-  - l)~- 


D R 
s 


where 


Ae . = e . 1 - e . 

i i i 


[19] 

[20] 


When  e,  and  e2  are  both  zero,  the  l/R  term  becomes  the  usual  coulomb  repulsion 
e2lC|,/D  R,  the  I/a,  term  the  well-known  Born  charging  term  for  ion  I 


(e('2/2a|)[(  I/D  )-l],  and  the  l/a2  term  the  Born  charging  term  for  ion  2. 


Step  2 


The  charges  are  now  given  by  Eq.  [21],  where,  v goes  from  0 to  I . 


e . ’ = e . ' + v(  e . -e . 1 ) 

i i ' ! \ ' 


[21) 


I if  i'j(r)  denotes  the  potential  at  the  end  of  step  I and  , ( r)  the  poten- 
tial at  any  State  of  step  2.  the  change  of  potential  during  step  ?.  for 
any  \j,  is  *'  ( r )->',(  r ) . Since  the  medium  responds  to  a change  of  charge 


J 


1 

i 


1 


J>e.v  only  via  the  optical  dielectric  constant  D , now,  we  have  during 


op 


step  2 


. Vi  v *>e  | 

'V  U = iw 


D r 


L22] 


Writing  5e  . ' as  v(  e . -e  . 1 ) and  V.  as  »■  J-Vj  we  have 


\n  = *i(rj  + 


-;(  e ,-e  , ' ) v[  e2-e,  1 ) 


D r , 
op 


D r 2 
op 


[23] 


where 


■Mr; 


+ 0 2 


Dr,  D r 2 
s s 


[24] 


$i  , the  average  potential  on  the  surface  of  ion  I minus  the  self-poten- 


tial is  obtained  by  subtracting  e|'/r,  from  Eq . [24]  and  v( e , -e  , ' j / r ( 

from  the  last  term  but  one  in  Eq.  [23],  then  replacing  r(  in  those  equa- 
tions by  a:  and  averaging  the  I / r 2 in  those  equations  over  the  surface  of 
ion  1,  thereby  replacing  l/r2  by  I / R Thus 


$ 


C t ^ I . U j 

a,  ' D 1 'DR 
s s 


v(  e i ~e  ■ ' ) ( J | , 

a , ' D 

op 


A e -e  1 i 

D R 
op 


[25 


The  work  done  during  this  step  is 


i 


W 


II 


- v . v r 

* i d e | t .'2  uc  2 

C v 0 


Vde, 


126] 


One  obtains 
W 


II 


ie^e,  e be  ,-e  , 1 Ae,  I ,1 

^ ^ -)-  + [r(Ae,)*(  — - Ij 


D 


op 


op 


6 1 5 a" 


[^1  Ae2)  2(  r-  I ) - e2'£e2(-!--  I ) ] --  . 


[27] 


op 


The  net  work  done  AG  is  the  sum  of  Wj  and  Wjj  and  is  the  free  energy  of 


this  fluctuation.  It  is  equa  I to 


WI  - WII 


'e? 


2a 


2 3 p 


n D D ' 
op  s 


[2  8] 


where  Ae . is  given  by  Eq . [20], 


If  AG^°'  denotes  the  free  energy  change  for  unit  concentration  of 


actants  at  a separation  distance  R,  it  is  related  to  the  same  quantity  at 
infinite  separation  AG0'  by  an  equation  similar  to  Eq . [I2l,  namely 


A_Gr 


A.G 


0.  . ,P  r 

r w - w 


On  the  intersection  hypersurface  of  Fig.  I,  the  leactants  and  products 
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have  the  same  distribution  of  configurations  (the  same  set  values  of  q*),  and 
the  same  potential  energy  (cf  Eq . [7]),  averaged  over  a distribution  of 
such  configurations.  Since  the  distribution  of  configurations  (and  mo- 
menta) are  the  same  for  the  reactants  and  products  on  the  intersection 
hypersurface,  the  entropy  is  also  the  same,  and  so  the  free  energy  of  the 
reactants  is  the  same  as  that  for  the  products  on  the  intersection  hyper- 
surface. It  then  follows  that  we  can  write 

AGr  - AGP  = AGft01  , [30] 

where  AG*"  is  given  by  Eq . [28]  and  AGP  by  the  same  expression  with  e.  re- 
placed by  e.P. 

To  find  e,‘  and  e2 1 one  minimizes  the  AGr  in  Eq.  [28]  subject  to  the 
constraint  imposed  by  Eq.  [30]. 

MGr  = (dAGr/^e  , ')  6e  , 1 -r  (oAGr/3Ae2 ')  5e2 ' =0  [31] 

6AGr  - 6AGP  = 0 . [32] 

Multiplying  the  second  equation  by  a Lagrange  multiplier  m and  adding  to 
Eq.  [31],  introducing  expressions  such  as  Eq.  [31]  for  5AGr  and  f^GP  into 
Eq.  [32],  and  setting  the  coefficients  of  5e , 1 and  be2 ' equal  to  zero, 
one  finds 


= e,  + m(e.-e.p) 


( 5 = 1-2)  . 


Introducing  this  result  into  Eqs.  [28]  and  [30]  one  obtains 


AG  — m \o 


where 


-(  2m+l ) \0  = AGR°‘ 


I I 


\o  ( Ae)  ( 2a  | * 2az  o 


and  Ae  is  e,p-elt  the  charge  transferred.  m is  the  solution  of  Eq . [35], 
The  free  energy  barrier  AG*  to  reaction  consists  of  the  work  term 
wr  to  bring  the  reactants  together  and  AGr.  Further,  solving  Eq.  [35]  foi 
m and  introducing  into  Eq . [34]  one  obtains 


AG*  w v (\o/4)[l  + UGr0,/\0)]  2 


In  the  electrode  case  the  electrostatic  potential  in  step  I is  given 


by  Eq.  [14],  but  with  the  e2  replaced  by  the  image  charge  of  reactant  I 
in  the  electrode  -e,.  The  image  charge  ensures  that  the  potential  given 
by  Eq.  [ 1 4 ] is  constant  on  the  surface  of  the  electrode,  where  r,  = r?. 
There  is,  in  addition,  another  term  due  to  the  interaction  of  ion  I with 
the  other  charges  on  the  electrode  and  with  the  surrounding  electrolyte. 
This  term  is  found  to  cancel,  apart  from  a minor  term,  in  steps  I and  2 
and  so  will  be  omitted  here  for  brevity.  The  derivation  continues  then 
as  before,  where  R now  denotes  the  distance  from  ion  I to  its  image, 
namely  twice  the  distance  to  the  electrode.  In  computing  Wj  and  Wjj  one 
computes  only  the  work  to  charge  ion  1.  One  finds  ultimately  expressions 
similar  to  Eqs.  [34] -[37]  but  now  we  have 

-(  2m+l ) \q  = - nF(  E — E 0 1 ) + wP  - wr  [38] 

and 

~ 2^  a7  ” D ' ^39' 

op  s 

That  is,  one  obtains  for  this  electrode  case 

AG*  = wr  + Uo/4)  [I  + [ -nF(  E-E 0 1 ) + wP  - wr}/\0]  2 [40] 

where  F is  the  Faraday. 

This  type  of  simplified  derivation  of  the  above  results  was  given 
earlier  in  Ref.  [12],  and  has  recently  been  made  more  readily  available 
i n a rev i ew  [13]. 

If  now  one  wished  to  include  s imu I taneous I y these  fluctuations  in 

solvent  polarization  and  those  described  earlier  in  reactants'  vibrations. 

One  could  readily  do  so;  To  the  right  side  of  Eq.  [28]  would  be  added 

Y £ k.(  q . -q  . °r) 2 , and  ~ T k.(q.-qj0p)2  wou  Id  be  added  to  the  correspond- 
i i 

ing  expression  for  AG  , One  would  then  proceed  as  before,  using  Eqs. 

[31] -[32],  but  now  Eq . [31]  would  contain  an  extra  term  TJ  k.(  q . -q  j r,V  ■ ‘q  - , 

i 

and  5AGP  would  contain  an  analogous  extra  term.  One  would  find  in  addi- 
tion to  Eq.  [33]  the  result  that 


i or 

q : = q : 


m(  q . °r -q  . oP) 


When  this  q.  and  the  e.1,  again  given  by  Eq.  [33],  are  introduced  into 
the  expression  for  AGr,  Eq,  [34]  would  again  follow,  but  with  X0  replaced 
by  \oxX|.  Xo  being  given  by  Eq . [36]  and  \.  by  Eq.  j 1 3]  . Similar  remarks 
apply  to  the  electrode  reaction  case. 
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DEDUCTIONS  AND  APPLICATIONS  OF  THE  THEORY 

When  a statistical  mechanical  treatment  is  undertaken  instead  of  a 
dielectric  continuum  one  again  obtains  equations  formally  similar  to  Eqs. 
[10]  to  [12]  and  [34]  to  [ 4 1 ] , but  with  Xo  having  a statistical  mechani- 
cal value  and  Xj+Xo  replacing  X-(  and  Xo  [3]  : 

= h Z exp(  -AG-*/kT)  , [42] 

where  in  the  homogeneous  case 

AG*  = wf  + (X/4)  [1+(AGR0,/\)]Z  - [43] 

AG0°'  is  given  by  Eq.  [29]  and  in  the  electrode  case  AGC°'  is  replaced  by 
K K 

AG°'  - -nF( E-E° 1 i + wp  - wr  , [44] 

K 

and  in  both  cases  \ is  tne  sum  of  two  contributions: 

X = Xj  + Xo  • [45] 

In  these  equations  k,  Z,  w1  and  wp  have  the  same  significance  as  before. 
AG°‘  is  the  free  energy  of  reaction  for  Eq . [4],  when  the  reactants  and 

products  are  each  in  unit  concentration  in  the  prevailing  medium  and  at 
the  prevailing  temperature.  E°‘  is  the  "standard"  half-cell  potential 
for  Eq.  [3]  under  the  same  conditions  and  E is  the  actual  half-cell  po- 
tential (in  tne  absence  of  concentration  polarization).  \ for  the  elec- 
trode case  differs  from  that  for  the  homogeneous  case,  as  before,  but  in 
each  case  contains  a vibrational  contribution  from  the  inner  coord i- 

i 

nation  shell  of  the  reactant  Xj  and  a contribution  Xo  from  the  solvent 
outside,  as  in  Eq.  [45].  Xj  is  of  the  form  in  Eq.  [13]  and  a dielectric 
continuum  estimate  of  Xo  is  of  the  form  in  Eq.  [36]  or  Eq.  [39], 

In  the  lecture  on  which  the  present  article  is  based  deductions  frot 
these  various  equations  were  described,  together  with  experimental  data 
regarding  them.  Since  similar  material  was  recently  presented  elsewhere 
[ 14]  , it  will  not  be  reproduced  here.  Instead,  some  of  the  deductions 
from  the  theoretical  equations  are  summarized  and  reference  is  rwide  to 
[14]  for  further  details  and  for  most  of  the  experimental  references.  Tie 
references  both  here  and  in  [14]  are  intended  merely  to  he  representa- 
tive. rather  than  complete.  Seveial  Ones  in  addition  to  [ 1 4 ] are  also  in- 
c 1 uded. 
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(1)  Cross-Reactions  and  Electron-Exchange  Reactions 

The  following  reactions  are  known  as  electron  exchange  reactions  and 
their  rates  are  most  frequently  measured  by  use  of  isotopic  tracers. 

A|(ox;  +A,(red)  - A |(  red;  + A,(ox)  [4o] 

A2(ox)  + A2( red)  A2(  red)  + A2(ox)  [47] 

Let  their  \'s  be  written  as  \M  and  X.22,  and  their  rate  constants  as 

k|t  and  k22.  The  reaction  in  Eq.  [4]  is  designated  as  a cross-reac- 

tion with  \ written  as  \|2  and  rate  constant  as  k,2.  Because  of  an 
additivity  property  [3]  of 

X 1 2 = ^ X ii  4 \2  z)  [4  8] 

and  the  form  of  Eq.  [43],  one  obtains  a prediction  of  k]2  in  terms 
of  k | | , k2  2 and  Km2,  the  equilibrium  constant  of  Eq.  [4],  when  the 
k's  for  the  individual  reactions  are  either  near  unity  or  nearly 
cancel,  and  when  the  work  terms  for  the  individual  reactions  are 
either  smail  or  nearly  cancel: 

k|2—  (k,tk22K|2f  2;  [49] 

where 

f 1 2 - (-f-n  K,2)*/4  k,  |k22/Z2)  . [50] 

Sctin,  referenced  in  [14],  pioneered  the  exper i menta I study  of  these 
relations  between  cross-reactions  and  se ! f -exchange  reactions. 

(2)  Dependence  of  the  slope  a of  a plot  of  in  kl2  vs  Tn  K|2  ( Bronsted 
coefficient)  (cf.  [14]). 

a = ^ ' + AGR°'/\)  . [51] 

(3)  Dependence  of  the  slope  j of  a plot  of  £n  i,  where  i is  the  current 
density  for  the  forward  reaction  in  Eq . [3],  s -nF(E-E0')  (Tafel 

s lope)  [ 14] . 

C ~[  I + [ -nF(  E-E 0 ' ) ^ wP  - wrV\]  [52] 


(4) 

(5) 

(6) 


Relations  between  electrode  rates  and  homoiieneous  reaction  rates 

[14]. 

Dependence  of  in  k^  on  solvent  dielectric  conti  rruurn  properties  when 
no  specific  reactant -sol  vent  interaction  exists  [ I 5] . 

Dependence  of  £n  k on  changes  in  equilibrium  bond  lengths  and  a'-;l«" 

[14]. 
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1(7)  Dependence  of  £n  k^  on  £e  (16). 

(8)  Formation  of  electronically-excited  states,  and  hence  chetni  lumines- 
cence, in  highly  exothermic  reactions,  and  properties  thereof  ( 1 4 ] . 

(9)  Contributions  to  the  entropy  of  activation  of  reactions  117]. 

(10)  Electrolyte  affects  in  electrode  reactions  [18]  and  in  homogeneous 
react  ions  [17], 

More  recently,  evidence  regarding  the  dependence  of  Tatel  slope  in 
electrode  reactions  on  E — E 0 1 may  be  found  in  [19]  and  on  chemi lumi nes- 
i cence  in  [20],  (For  further  theoretical  studies  or.  references  related 

] to  highly  exothermic  reactions,  such  as  those  involving  chemiluminescence, 

see  [14].)  A nice  example  of  a test  of  Eq.  [43]  is  given  in  [21],  An 
application  of  the  present  electron  transfer  theory  to  electron  transfer 
reactions  with  negative  activation  energy  [22]  is  found  in  [23],  Other 
interesting  developments  in  electron  transfer  theory  include  the  study  of 
reactions  at  semiconductor  electrodes  [24-26], 

COMMENTS  ON  SEVERAL  DEVELOPMENTS  SINCE  1957 

Elsewhere  in  this  volume  OMR  Technical  Report  No.  12  ( 1957)  is  re- 
produced. It  contains  the  formulation  of  this  electron  transfer  theory 
for  electrode  reactions.  That  Report  was  concerned  with  the  dielectric 
continuum  contribution  to  the  free  energy  barrier,  a contribution  for 
which  a simplified  derivation  was  given  in  an  earlier  section  of  the  pre- 
sent paper.  The  derivation  in  ONR  Technical  Report  No.  12  is  given  in 
terms  of  vector-fields,  vectors  associated  with  the  dielectric  polariza- 
tion, with  the  electric  field  due  to  the  charges  themselves,  and  with  the 
total  electric  field  due  to  the  charges  and  the  polarization.  In  the 
simplified  derivation  given  in  a previous  section  the  electrostatic  quan- 
tities are  not  expressed  in  terms  of  vector-fields  but  rather  in  terms  of 
scalar-fields  associated  with  potentials  and  charges,  though  fields  of  a 
restricted  kind,  namely  of  the  functional  form  given  by  Eq.  [I4j.  A more 
general  derivation  is  given  in  terms  of  sea  I ar-f i e I ds  in  [3],  and  vector- 
fields  in  ONR  Technical  Report  No.  12.  No  assumption  is  made  initially, 
for  example,  that  the  reactant(s;  is(are)  spherical. 

One  genera  I izat i on  of  the  dielectric  continuum  result  in  Eqs.  [36] 
and  [39]  for  \0  is  to  include  the  fact  that  there  is  a dispersion  of  fre- 
quencies characterizing  the  dielectric  response  |27,28],  The  correct io> 
in  the  \0  value  amounts  to  18  [28). 

One  point  Of  interest  regarding  derivation  in  [9,10]  terms  ot  vectoi - 


fields  (viz.  electric  vectors)  versus  one  in  terms  of  scalar  fields  (viz. 
charge  densities  and  electrostatic  potentials)  in  [3]  is  the  following: 

In  the  vector-field  derivation  in  QNR  Technical  Report  No.  12  [10]  and  in 
its  predecessor  [9]  the  transcription  of  the  final  vector-field  equations 
for  the  electrostatic  free  energy  into  scalar-field  quantities  was  needed 
and  was  made.  This  transcription  is  avoided  by  using  instead  scalar- 
field  expressions  throughout. 

Two  of  the  generalizations  in  Ref.  [3]  are  to  use  statistical  mech- 
anics instead  of  dielectric  continuum  theory  and  to  treat  the  inner  shell 
contribution  to  the  fluctuations  in  coordinates.  ^This  second  generali- 
zation was  given  also  in  Ref.  [29].)  The  generalizations  led  to  similar 
predictions  as  those  given  earlier  in  the  present  paper,  but  without  the 
simple  dielectric  equations  [36]  and  [39]  for  X0  and  hence  without  the 
simple  predictions  of  solvent  effects  for  reactants  which  do  not  specifi- 
cally interact  with  the  solvent.  A statistical  mechanical  expression  for 
was  given  instead, [3],  one  whose  evaluation  must  await  further  appli- 
cation of  the  statistical  mechanics  of  polar  liquids.  When  a dielectric 
continuum  estimate  was  made  of  this  contribution  one  obtains  the  result 
for  \o  given  in  Eq.  [36]  for  homogeneous  reactions  and  Eq . [39]  for  elec- 
trode reactions.  When  the  inner  contribution  is  included,  the  same  final 
equations  obtained,  but  with  \0  replaced  by  \0+Xj,  as  in  Eq.  [45]. 

In  the  statistical  mechanical  treatment  it  was  recognized  [3,29]  that 
when  the  q in  Fig.  1 in  the  vicinity  of  q*  is  primarily  a solvent  orienta- 
tional or  solvent  vibrational  coordinate  the  free  energy  change  giver  y 
Eq.  [43]  is  actually  the  free  energy  of  a fluctuation  of  coordinates  to 
values  centered  on  the  intersection  hypersurface  rather  than  being  con- 
fined to  it.  This  effect  introduces  a minor  correction  factor  o [3]  in 
the  pre-exponential  factor  in  Eq.  [42] ; p is  of  the  order  of  unity. 

Among  other  developments  since  1957  has  been  the  recognition  ot  t In- 
various  consequences  of  the  theoretical  equations,  listed  earlier.  Tit 
relation  given  by  Eq.  [49]  was  first  given  essentially  in  I9o0  [29],  ap- 
plied to  experimental  data  in  1963  [30],  and  given  a more  general  deri- 
vation in  1965  [3].  The  dependence  of  Tafel  coefficient  on  E-L  ' was 
not  ronf  i r rried  until  197  5 [19],  and  use  was  niade  of  Mohilrei's  treat- 
ment [18]  of  electrolyte  effects  to  calculate  the  work  terms  w'  ard  w1 
The  solvent  effect  was  tested  in  1970  [15]  and  the  effect  ol  £e  ii  Eq 
[39]  in  1969  [ 1 6] . A comparison  between  some  rates  of  electrode  react i o>  s 
and  the  electronic  structure  of  the  reactin'  species,  aimed  at  correlat- 
ing qualitatively  with  charges  of  equilibrium  I ond  lengths  is  given  in 
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Ref.  [14].  Applications  of  the  cross -re  I at i on  Eq.  [49]  to  biological 
systems  have  also  been  made  and  referenced  in  [14], 

Eq . [49]  has  been  widely  tested  (referenced  in  [14])  and  has  been 

very  useful  in  correlating  a large  body  of  data.  The  main  anomaly  ob- 
served thus  far  in  this  equation  is  for  reactions  of  Co+3(aq)  ion.  This 
ion  undergoes  an  extensive  electronic  rearrangement  in  forming  Co+2(aq) 
and  a possible  explanation  for  the  anomaly  is  given  in  [14], 

Quantum  effects  on  the  vibrational  or  solvent  motion  have  also  been 
treated  [7,8,13,31b],  At  sufficiently  low  temperatures  the  system  does 
not  use  the  energy  AU*  in  Fig.  I to  go  from  the  R curve  to  the  p curve. 
Instead  it  can  "tunnel"  through  the  barrier  there.  This  tunneling,  in 
many-d i mens iona I q-space,  has  been  treated  by  means  of  Fermi's  Golden 
Rule  for  rad iat i on  1 ess  transitions  and  Franck-Condon  factors  [ 7, 8, I 3, 3 1 b] . 
At  high  temperatures  the  expression  for  the  free  energy  barrier  reduces 
to  the  classical  one  given  earlier  in  this  paper.  An  analogous  formalism 
has  been  used  to  treat  [31]  biological  electron  transfers  [32]  at  very 
low  temperatures. 

In  the  quantum  treatment  [7,8,13]  a treatment  of  the  violations  of 
the  reactants  is  relatively  straightforward,  but  the  solvent  vibrations 
are  treated  as  though  they  are  vibrations  of  a so!  in  and  hence  undamped. 

In  the  liquid  there  is  a strong  damping.  (This  approximation  was  not 
made  in  the  classical  statistical  mechanical  treatment  [3i.)  A modifi- 
cation for  the  damping  on  the  quantum  mechanical  result  nas  been  des- 
c r i bed  [33] . 

Again,  when  curve  P in  Fig.  I crosses  curve  R on  the  left  hand  side, 
so  that  the  slopes  of  both  curves  are  both  negative  at  the  crossing  point, 
the  transition  from  R to  P curve  can  only  occur  nonadi abat i ca II y and  quan- 
tum treatments  involving  the  Franck-Condon  factors  and  related  closely  to 
the  just  mentioned  quantum  formalism  [7,8.13],  have  been  used,  referenced 
in  [14].  Such  treatment  are  needed  for  highly  exothermic  reactions'  wmi" 
AU  in  Fig.  i is  sufficiently  negative  the  crossing  will  indeed  be  such 
that  the  R and  P curves  have  the  same  slopes  at  the  crossing  point. 

Finally,  there  are  reactions  of  the  redox  or  electrochemical  type 
where  bonds  ore  actually  broken  and  formed  during  tire  reaction.  Here, 
the  description  of  the  potential  energy  surfaces  given  in  Fig.  I is  sug- 
gfestive  though  not  adequate.  Other  models  have  been  used,  the  bond 
energy-bond  order  model  [34]  for  example  [3b, 36],  and  an  analog  to  Ecs. 
[43],  [40]  and  [!>0]  has  been  obtained  [' 


but  will. 


where 


AG*  = wr  + ( X/4)  + (AGR0'/2)  + ( AGR°'/2y)/.n  cosh  y 

f j c°sh  y)/y 

'12  - t I 2 > 

y = (2VAGr°'/\)  = (In  Kl2)Y/Tn(k,  ,k22/Z2)  , y = Tn  2 


(The  case  of  the  hydrogen  discharge  reaction,  Eq.  [I],  will  oe  treated  in 
a forthcoming  publication  [I]  using  a combination  of  nonequilibrium  po- 
larization, BEBO,  and  Franck-Condon  overlap  methods  to  treat  different 
aspects  of  the  overall  problem.)  Whereas  the  slope  of  a Bronsted  plot 
(homogeneous  reaction)  or  Tafel  plot  (electrode  reaction)  was  given  by 
Eqs.  [ 5 1 ] and  [52],  they  are  given  [36]  in  a BEBO  model  by  Eq . [56]  when 

the  work  terms  can  be  neglected. 

a = |(  I f tanh  y)  1 56] 

where  y is  given  by  Eq.  [55].  In  the  electrode  case  the  2AG  °'/\  is  re- 
placed  by  2{-nF(E-E0')  + wP  - wr]/X.. 

The  theory  of  electron  transfer  reactions  in  solution  and  at  elec- 
trodes continues  to  be  in  an  active  and  developing  state.  An  estimate 
has  been  made  for  the  h in  Eq.  [42]  for  the  fer rous - fer r i c electron  ex- 
change reaction  [37].  Calculations  are  desirable,  using  increasingly  ac- 
curate theories  of  electronic  structure,  for  x's  of  this  and  other  reac- 
tions including  those  of  Co43(aq)  for  which  extensive  electronic  rear- 
rangement occurs.  Such  calculations  will  also  permit  more  insight  into 
orientational  effects  [38]  in  electron  transfer  reactions.  Again,  dev- 
elopments can  be  expected  in  statistical  mechanical  evaluation  of  the 
term  obtained  for  fluctuations  in  orientations  of  solvent  molecules, 
needed  to  reach  the  intersection  region  of  Fig.  I.  There  have  been 
numerous  enlightening  interactions  between  theory  and  experiment,  and  we 
can  continue  to  look  forward  to  this  fruitful  interplay  in  the  future. 
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